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Abstract. The critical exponent corresponding to the renormalization of 
the composite operator i/jip is computed in quantum electrodynamics at 
0(1/Nj) in arbitrary dimensions and covariant gauge at the non-trivial zero 
of the /3-function in the large Nf expansion and the exponent corresponding 
to the anomalous dimension of the electron mass which is a gauge indepen- 
dent object is deduced. Expanding in powers of e = 2 — d/2 we check it is 
consistent with the known three loop perturbative structure and determine 
the subsequent coefficients in the coupling constant expansion to all orders 
in MS. 
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The fundamental functions of the renormalization group equation, such 
as the /3-function and mass anomalous dimension, are crucial ingredients 
in determining the properties of the Green's functions of a renormalizable 
quantum field theory and how they depend on the energy scale. They are 
ordinarily determined by renormalizing Green's functions at successive or- 
ders in perturbation theory where the coupling constant, g, is assumed to be 
small. Invariably one finds that calculations become exceedingly difficult at 
a certain order, due to the large number of Feynman graphs to be evaluated, 
which is therefore tedious, and the intricate nature of some of the integrals 
which can arise. Thus the renormalization group functions are only known 
to the first few orders for most theories. For example, the /3-function of QED 
has recently been calculated to fourth order in MS in [|IJ by use of computer 
algebra packages and the wave function and mass renormalization functions 
are known to third order, by restricting the QCD results of 0] to the abelian 
case, which extended the earlier two loop calculations of |3| Clearly it is 
important to have methods which circumvent these calculational difficulties 
and allows one to gain an insight into the perturbative structure of these 
functions to very high orders in the coupling constant. 

One way to achieve this is by computing the functions in terms of a 
different expansion parameter. For example, in theories which possess N 
fundamental fields, such as Nf electrons in QED, one can use 1/N as an alter- 
native expansion parameter. In the earlier work of || and Q the /3-function 
and mass renormalization, which are gauge independent quantities jj], S], 
were computed at leading order in large Nf in QED. This was achieved by 
explicitly computing the simple poles of the leading order Green's functions, 
in dimensional regularization, where the photon propagator was replaced by 
a chain of one loop bubble graphs. The simple poles were then absorbed 
minimally into a renormalization constant and the (4 — 2e)-dimensional per- 
turbation series of f3{g) and r y rn {g) were deduced, which gave the 0(1/Nj) 
coefficients of each function to all orders in g. Whilst this was successful in 
giving agreement with what had previously been calculated in MS and also 
provided new coefficients, it is certainly not easy to extend that analysis to 
probe the perturbation series more deeply. 

An alternative approach is to use the large Nf self consistency method 



introduced in [g, |10[ where d-dimensional critical exponents are computed at 



the non-trivial zero of the /3-function. As the theory is finite and possesses a 
conformal symmetry there, the renormalization group equation simplifies to 
the extent that the critical exponents one computes are related to the critical 
renormalization group functions which also allows one to extract all orders 
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information. The beauty of applying the method to a four dimensional gauge 
theory such as we do here for QED is that by using dressed propagators one 
substantially reduces the number of Feynman diagrams to be analysed and 
crucially it is possible to achieve results at 0(1/Nj). Previously using these 
techniques in QED the electron wave function exponent r\ was deduced at 
0(1/ N f ) in the Landau gauge in Q and later at °( l l N f) in El- The 



original /^-function calculation of || was reproduced more concisely in 
at 0(1/Nf) and the gauge independence was checked in [14] by computing in 
an arbitrary covariant gauge. The significance of the agreement of [13] with 



is that it establishes that the exponents calculated by the self consistency 



method of ||, 10] encode MS information. 



In this letter, we extend the work of || by computing the electron mass 
anomalous dimension, "y m (g), to 0(1/N 2 ) by determining the relevant crit- 
ical exponent. Its calculation is necessary as an independent check on the 
three loop result of || as well as providing the remainder of the perturba- 
tion series for j m (g) at 0(1/Nj). This is required for future four loop QED 
and QCD calculations and will also serve as the foundation for finding the 
0(1/ Nj) quark mass anomalous dimensions. Simultaneously we deduce a 
result for three dimensional QED which is important since it will provide 
improved estimates for a gauge independent quantity which can be measured 
on the lattice. This will give a useful check for the Monte Carlo algorithms 
used to simulate problems such as the existence or otherwise of a chirally 
symmetric phase. Finally, we draw attention to the fact that the procedure 
we follow here is based very much on the 0(1/ N 2 ) calculation of the mass 
operator in the 0(N) Gross Neveu model, fl5| . 

To compute critical exponents in large Nt in QED we use the lagrangian 



in the form, [11], 



4e 2 2£e 2 



#w + A^yy -t%- -^t^- w 



where V* is the electron field, 1 < i < Nf, F^ IU = d^A v — d u A^, £ is the 
covariant gauge parameter and is the U(l) gauge field. The electron 
coupling constant, e, has been rescaled into the gauge field kinetic term to 
ensure that the vertex has the correct structure for applying the method of 



uniqueness of [16] which is the main technique for computing the massless 
Feynman integrals which arise. The location of the critical coupling, g c , in 
whose neighbourhood we will solve (1) is defined as the non-trivial zero of 
the d-dimensional /3-function. It has been calculated in MS to fourth order 
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in pi in (4 — 2e)-dimensions using dimensional regularization as 



0(g) 



(d - % 



2N f 2 N f 3 

~t 9 + i 9 



N f (22N f + 9)_ 4 



64 



616jV/ / 416C(3) _ 380 \ 
243 + V 9 ~27~) 



144 

N f + 23 



g 5 + 0(g 6 ) (2) 



where we use the conventions of |l7j but have defined g = (e/2ir) 2 . Thus 
to deduce information on the (4 — 2e)-dimensional perturbation series from 
the exponent we will calculate, (2) implies that then 



9c 



3e 



27e 2 



+ 



99e d 



+ 



77e 4 



AN 2 IQNf IQNf 



(3) 



In the critical region the field theory has a conformal symmetry and so the 
propagators and Green's functions obey a simple power law behaviour. For 
instance, in momentum space, as k 2 — > oo, the propagators of (1) satisfy, 
0, 



fi—a 



Ai/a (k) 



B 



(k 2 y~p 



i 



. k u k, 



k 2 



(4) 



where A and B are the respective ^-independent amplitudes and a and (3 
are the critical exponents of each field and their canonical and anomalous 
dimensions are defined to be, W%, 



a 



1 + iT? 







i-v-x 



(5) 



where rj is the electron anomalous dimension and x is the anomalous dimen- 
sion of the 3- vertex. Each is 0(1/Nf) and depends on d = 2/x. The former 
corresponds through an examination of the critical point RGE to the elec- 
tron wave function renormalization which is known to third order in and 
r\ has been deduced at 0(1/N 2 ) in the Landau gauge in [12]. Although it 
is a gauge dependent quantity it must always be calculated first within the 
self consistency formalism of || as it plays a central role in the computation 
of other exponents. 

To compute the electron mass exponent, 'y m (g c ), we follow the same 
strategy as the two and three loop perturbative approaches of [2-4]. There 
7m(<50 w as calculated indirectly by first renormalizing the composite operator 
"0*0 and then using 

7m(g) = 1^(9) + 72(5,6 (6) 
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where 72(9, £) is the wave function renormalization constant which differs in 
overall sign to the definition given in ]|]. In exponent language (6) translates 
into 

7m(Sc) = ?7 + 7^(ffc) (7) 

Since both r\ and 7^(g c ) are gauge dependent the cancellation of the pa- 
rameter £ when we calculate with an arbitrary covariant gauge will provide 
a very stringent check on our computation before checking with [||. For 
completeness, in the same notation as (2), M, 



7m (fiO 



3 (20iV> - 
^ + (46- 



%2 



48C(3))iV) 



129 



32 



+ 0(g 4 



(8) 



We will now detail the leading order analysis to reproduce the result of 
H but in the critical point self consistency approach which will also illus- 
trate the simplicity of the method and will allow us to elaborate on several 
technical points in preparation for the 0(1/Nj) calculation. First, we have 
to introduce a regularization to handle the infinities which will arise in the 
calculation. This is achieved by shifting the exponent x by X ~^ X + A where 



A is an infinitesimal quantity, [12]. Once this is introduced the Green's func- 
tions which one calculates in (^-dimensions will have a three part structure, 
Ipm , |Ts| l. It will involve terms which have poles in A which are removed 
by a conventional renormalization, |j[q| . This leaves O(l) terms which are 
multiplied by hip 2 , where p is the external momentum, and O(A) terms 
which vanish as the regularization is lifted. As we are working at a critical 
point where p 2 — > 00 the presence of hip 2 terms will spoil the scaling be- 
haviour. General arguments, however, show that their sum will be of the 



form (p 2 )2 7 V"/'( 5lc ) ) [18, 19]. Therefore, expanding this in powers of 1/Nf and 



choosing 7<^(g c ) appropriately at each order will remove these hip terms 
to leave a Green's functions with sensible scaling behaviour. Indeed this 
observation which was discussed originally in [|l§|] was also used to renor- 



malize the 2-point function of QED in [12, 2C[. In other words to deduce 
l^ipidc) a t each order in 1/Nf one computes the Feynman diagrams with 
dressed propagators which contribute to the inclusion of ipi/j as a composite 
operator in some Green's function, which for this letter will be {ip[4>vp]'ip) , 
isolates the Inp 2 terms after A-renormalization and sums their coefficients. 

At leading order there is only one Feynman graph to consider which is 
given in fig. 1. The circle with a cross indicates a xjjijj insertion. Substituting 
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the asymptotic scaling forms of the propagators of (4) for the lines of fig. 1 
and computing the trivial one loop integral, it is 

z(2n - 1 + £)a(ji - 1 + A)a(l)a(/i - A)(p 2 )- A (9) 

where we have defined a(a) = r(/j, — a)/T(a) for all a and set z = A 2 B. 
Expanding (9) in powers of A one obtains 



(2 M - 1 + j) Zl 

Ar(/x)iV) 



l-Alnp / + 7 +0(A') 



(10) 



where z = Yh^Li Z i/Nj and z\ is deduced from the leading order computation 
of T) in [llj] as 

Z1 = ?m (in 

Thus, 

/z(2/x-l + £K 
7^,i(ffc) - - ^32)(2^i) ^) 

o (2//-l)fc-2)r(2//) 

7,1 4r2(^)r( / x + i)r(2- / ,) 1 J 

which corresponds to the Landau gauge expression for r/i. It is a trivial 
exercise to repeat the computation of 771 for non-zero £ to deduce 

* _ [(2/, -!)(/, -2)+^° 
771 (2^-1)0* -2) 1 ' 



where 



which from (7) implies 



/ x 2 »7l 

lmA9c) = ~ y^y (15) 



in agreement with Q, but deduced by considering only one Feynman graph. 

The way to proceed to 0{1/Nf) is now clear and involves first of all 
computing the contributions from fig. 1 since its lines are represented by 
propagators which possess Nf dependent exponents as well as the obvious 
counterterm graphs. However, the main effort lies in considering the higher 
order graphs of fig. 2. The fourth graph can be determined by a simple 
extension of the one loop integral of fig. 1 and it contributes 

3^-1 + ow 



2(/i -!)(//- 2)2(2^-1)2 
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to T^vi 2(9c) ■ To calculate the remaining graphs one first maps each integral 
into coordinate space through the Fourier transform 



2\a 



Akx 



a(a) 



(17) 



Since the third graph of fig. 2 then corresponds to an insertion on the com- 
pletely internal electron line of the second order correction to the electron 
self energy it is A-finite in coordinate space, |pQ |. Also, the values of the first 
two graphs are equivalent and after mapping to coordinate space they cor- 
respond to a self energy graph where the insertion is now on an electron line 
joining to an external vertex. It is not instructive to reproduce the tedious 
algebra connected with the calculation of these graphs aside from remark- 
ing that the techniques of integration by parts, |l4|j , subtractions, [1C], and 
uniqueness, fl6[| , we used have already been discussed in the computation of 
the electron self energy at second order, [|l^, ppR . After a substantial amount 
of algebra we found that the first graph of fig. 2 contributes 



0\2 



(2^-1)^-2)! 



(2//-l + fl[(2/x-l)fc-2)+^] 

(m-i) 

(2^-l)( M -l)(l-£) 



(18) 



to (.9c) 1 whilst the third gave 



2( M -l)(2/x-l)2( M -2)2 



[(2^-1 + 0(2^-3-^-4^-1)] (19) 



The contributions from the graph of fig. 1 which involves the various coun- 
terterm pieces are 



M (2 M -l + OW) 2 



(2 M -1)(3a*-4) 



(/i-l)(2/i- 1)2(^-2)2 

+ 3£/i - 3n(p - l) 2 ( e 



(20) 



where 0(/i) = V' / (/ i ~~ 1) — ^'(1)) VK 3 -) is the logarithmic derivative of the 
T-function and we have used Xi = ~ "Hii [0]> an d 



22 



2(/i-2)2(2 M -l)* 



»2 r 



e 



1 



(21) 
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from the x-space calculation of 772, [12|. Thus, we have 



lifnf>,2 (dc 



0\2 



(2^-1)2(^-2)2 
2/x(2/x - 1) - 3/i(2/x 



2(1 -Q(2/i -!)(//-!) 
A* 

-i + e)( M -i)fe 



(m-i) 



(22) 



Finally, to obtain 7 m ,2(5c) we need to add 7/2 in an arbitrary gauge to 
(21). We have deduced this by adapting the Landau gauge computation of 
[2l to find 



'/2 



1)(^?) 2 



(^-2)2(2/^-1)2 
+ 3/i[(2/i-l)(/i- 



2(2/i - 1) 
A* 

2) + £m] 



[0* 
e - 



•1)0* 
1 



3) + ^] 



From (7), (22) and (23) then we have simply 
6(r/?) 2 



7m,2(ffe) 



(/i -2)2(2/1-1) 



/i(/i- 1)6 



(m-i) 



1 4/i 
+ - + — - 2 
/x 3 



(23) 



(24) 



where it is reassuring to note that, as in (15), all ^-dependence has cancelled, 
which is the first check on our result. 

The subsequent check is to set /1 = 2 — e in (24) and expand in powers 
of e to 0(e 3 ) in order to compare with (8) at the value of (3). We find 
that both (24) and (8) are in agreement and note that the £(3) of (8) arises 
from the expansion of in (24). Thus the gauge independent result (24) 
corresponds to the 0(1/Nj) correction to (15). Moreover, one can now 
deduce new coefficients which will appear in "f m (g) at higher orders in g by 
expanding ^ m ,2(9c) to the subsequent orders in e. The 0(1/ Nj) corrections 
for g c are required for this and they are encoded in the /^-function exponent 
A = 



Wi9c) of 



A 



(2/i -3) (/1 



f 



(25) 



in the present notation. We record, for example, that in MS to fourth order 
in g, (8) will become 

„3 



im(g) 
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'WN f 3" 
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-V 144 " C(3) + -A 5C(3) - 3C(4) - 04 +pNf + Q 



g 4 (26) 
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where the unknown coefficients p and q can only be determined by an ex- 
plicit 4-loop calculation or respectively knowledge of 0{1/Nj) and 0(1/Njr) 
exponents. However, it is worth pointing out that one can now make use of 
(26) to substantially reduce the amount of work one would have to do in the 
explicit calculation. By isolating only those Feynman graphs with one and 
no electron loops the residues of the simple poles in e respectively determine 
p and q. 

Finally, since (24) is valid in d-dimensions we can restrict to fx = 3/2, to 
obtain 

*.<*)- -^-^I^-^+O^) (27) 

As this also is a gauge independent quantity it can be used to determine 
exponents for relatively low values of Nf to compare with numerical work. 
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Figure Captions. 

Fig. 1. Leading order graph for 7^,(<7 C )- 
Fig. 2. Higher order graphs for 7^ (<?<;)• 
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